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Abstract. We consider quantum dynamical systems specified by a unitary 
operator U and an initial state vector if>. In each step tlie unitary is followed 
by a projective measurement checking whether the system has returned to 
the initial state. We call the system recurrent if this eventually happens with 
probability one. We show that recurrence is equivalent to the absence of an 
absolutely continuous part from the spectral measure of U with respect to (j>. 
We also show that in the recurrent case the expected first return time is an 
integer or infinite, for which we give a topological interpretation. A key role in 
our theory is played by the first arrival amplitudes, which turn out to be the 
(complex conjugated) Taylor coefficients of the Schur function of the spectral 
measure. On the one hand, this provides a direct dynamical interpretation of 
these coefficients; on the other hand it links our definition of first return times 
to a large body of mathematical literature. 

1. Introduction 

One of the standard questions addressed in the theory of Markov chains is the 
distinction of transient and recurrent processes. By definition, a Markov chain is 
called recurrent, if almost every path comes back to its starting point, otherwise 
transient. In this paper we discuss an analogous question for unitary quantum 
dynamics. Immediately, this presents us with a problem, which is typical for many 
generalizations of classical concepts to the quantum world: The definition as given 
for Markov chains clearly requires some monitoring of the process: we have to 
check after every step whether the particle has returned. But this monitoring, if 
it is to give any non-trivial information about the system, necessarily changes the 
dynamics. Therefore, there are two options: We can either try to reformulate the 
problem in such a way that the monitoring is not needed, or else we include the 
monitoring into the description. An approach of the first kind has recently been 
put forward by Stefahak, Jex and Kiss |33j on the basis of a formula, which Polya 
derived for Markov chains, and which is only based on the return probabilities, i.e., 
the probability for the particle to be back at a given time, without any conditions 
on the path leading there. We call unitaries with specified initial state satisfying 
the criterion 33J SJK-recurrent. In this paper wc choose the second option. Our 
notion of recurrence implies SJK-recurrence, but not conversely. We refrain from 
comparing our definition with a host of other ones in the literature. 
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A common feature of these approaches is that the recurrence criterion depends 
only on the unitary operator U and the initial state 0, and hence only on the scalar 
measure fj,{du) — {(j)\E{du)(j)) on the unit circle, which is obtained by the spec- 
tral measure E oi U (so U — J uE{du)). Ultimately, both criteria have a simple 
expression in terms of properties of /i. This connects the recurrence problem to a 
rich mathematical literature about measures on the unit circle and their orthogonal 
polynomials (the canonical review of the field is [301 IS])- The first return prob- 
abilities in our approach are the squared moduli of the Taylor coefficients of the 
so-called Schur function [2S] of the measure, which so far did not seem to have a 
direct dynamical interpretation. Our main result is that the process is recurrent iff 
the Schur function is "inner", i.e., has modulus one on the unit circle. Furthermore, 
we show that the winding number of this function has the direct interpretation as 
the expected time of first arrival, which is hence an integer. 

We came to this investigation through our respective work on quantum walks 
[H [S HI ■ Quantum walks which are basically translation invariant (up to a local 
perturbation) are not recurrent, but the analysis of first returns is by no means 
limited to deciding recurrence. We show how the methods developed in the quantum 
walk context can be used to get explicit expressions in simple cases. 

Our paper is organized as follows. In Sect. [2] we recapitulate the classical idea 
of recurrence and the quantum version by |33| . and then propose our alternative 
based on monitoring the dynamics. In Sect. [3] we give a spectral characterization of 
this notion (Theorem [ij and, for the sake of comparison, a similar characterization 
for the SJK-recurrence. The determination of the expected first return time is in 
Theorem [2] of Sect. |4] the main result is that the expectation of the first return time 
is, surprisingly, an integer or infinite. Again detailed results on the corresponding 
SJK-quantity are provided for comparison. Sect. [5] treats in more detail the finite 
dimensional case and shows, in particular, how the variance of the first return time 
diverges when this integer value changes. Examples from the field of quantum walks 
are presented in Sect. [6] We provide two rather different methods for computing 
return probabilities, which however, have a large overlap of applicability. We close 
with an outlook on some examples of singular spectrum and natural generalizations 
of our basic notion of recurrence. 

2. Notions of recurrence 

In this section we compare our notion of recurrence to the one introduced by 
Stefanak, Jex and Kiss |33j . which avoids introducing explicit monitoring steps. 
Their idea is based on the theory of Markov chains, particularly the idea of recur- 
rence as initiated by Polya [3S]. We briefly review this classical background [TS] 
here, partly to improve the understanding of [33^, but chiefly as a basis to point 
out close analogies between the classical return probabilities and return amplitudes 
arising later in our paper. After introducing our proposal we provide an example 
demonstrating that the two approaches may come to opposite conclusions. 

2.1. The classical case. A (discrete time) Markov chain on a countable state 
space X is defined by its transition probability matrix Pxy, where x,y & X. From 
the interpretation of Pxy as the probability for the process to move from x to y in 
one step it is clear that Pxy > and J^y Pxy = 1- Then the Markov property entails 
that the probability for going from x to y along intermediate steps (xi, . . . , a;„_i) is 
Px,xxPxi,x2 ■ ■ ■ Px„-i,y If we do not care for the intermediate steps these variables 
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are summed over, resulting in the matrix power P"j^. Now we single out a point 
€ X. It will be the initial point for all paths we consider, and we are interested 
in the returns to 0. The probability to return in n steps is obviously 

(1) Pn—Pm— ^ Po,XiPxt,X2- ■ ■ Px^^ifi- 

Xi,. ..^Xn-1 

We would also like to know the probability g„ for this same transition subject to 
the condition that no intermediate step is 0. That is, 

(2) qn ^ Po,xiPxt,X2 ' ' ' Px^-xfi 

Xi ....,Xn-l 

aU^o 

These numbers are the probability law for the random variable first return time. 
Clearly, ^ is the sum over just a subset of the terms in (jlj. If we split the sum 
([T]) according to the highest index k at which Xk = 0, we get a sum in which the 
summation indices with £ < k are unconstrained, but those with i > k must be 
non-zero. In other words 

n 

(3) Pn = y^^Pkqn-k, 

fc=0 

where for the endpoints we use the natural conventions po — 1 and qo = 0. This 
formula is not valid for n = (where it would give the wrong equation 1 — pq — 
Poqo — 0). Therefore, if we introduce the generating functions 

OO CO 

(4) p(z)^^p„z" and 9(2) = ^ g„z" 



n=0 Ti=0 



and sum over ([3| (the n = term has to be treated separately) we get 

(5) piz) = l+p{z)q{z). 

Hence the return probability or, more precisely, the probability to eventually return 
to the initial state 0, is now 

(6) i?^ = ^ g„ = g(l) = 1 - 4t - 1 - • 

Therefore the Markov chain is recurrent from the point if and only if i?*^ — 1, or 

2.2. Recurrence without monitoring. Now the criterion for recurrence pro- 
posed by j33j is to apply ^ to the sequence of quantum return probabilities 

(7) Pn = |(0|C/»P. 

So a unitary U with initial state (j) is called SJK-recurrent iff X^n^*" = 00. It is 
clear that the entire computation leading up to the criterion in the Markov case, 
being based on conditioning, loses its significance in the quantum case. Also the 
operational meaning as the probability of some property of paths of the process 
is lost. Therefore [S3] proposes an alternative operational criterion, based on the 
observation that J2nPn = 00 is equivalent to Yini^ ~ Pn) = 0. [33] defines the 
return probability i?^''^ by means of 

00 

(8) R^''' = l-l[{l-pn). 
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Operationally, this corresponds to the following procedure [33] : "Take a system and 
measure the position of the walker after one time step at the origin, then discard 
the system. Take a second, identically prepared system and let it evolve for two 
time steps, measure at the origin, then discard the system. Continue similarly for 
arbitrarily long evolution time. The probability that the walker is found at the 
origin in a single series of such measurement records is the Polya number [i.e., 
jjSJKj ii T]^is procedure is strange indeed. Imagine what it would mean in the 
parable of the prodigal son. The father is clearly interested in the return of his son, 
but when he does not come back until year n, he sends out another son and asks 
whether this one returns exactly in year n + 1, ignoring all earlier and later arrivals 
by this son, and indeed of all those sent before. While being rather wasteful of 
sons, it still does not provide an answer to the question "Does every prodigal son 
eventually return?" . Yet the procedure of |33i is probably the best we can do if 
we want to avoid observing the quantum process. We will avoid such strangeness 
in our approach but, of course, this means that we must disturb the free evolution 
and explicitly describe the measurements. 

2.3. Arrival by absorption. Here we go the alternative route of explicitly intro- 
ducing measurements to check the return of the system after every step. Starting at 
00 = ((>, one step thus consists of the unitary transformation 0„ U4>n followed by 
the measurement of the projection |(/))((/)|. If this gives a positive result (for which 
the probability is \{(j)\U(j)n)\'^) the experiment is over. Otherwise the system is left 
in the state </>„+i = c„+i(I — |0) {(f)\)U<pn, where we choose c„+i so that ||0n+i|| = 1. 
In this iteration we will always have 0„ proportional to where 

(9) U=il-\^){^\)U. 

Hence we can set (j>n = t/"<^/||?7"0||. The normalization factor s„ = ||J7"0|p has 
the direct interpretation as the survival probability, i.e., for the walker to remain 
undetected for all steps including the n^^. The probability for detection in step 
n, conditioned upon survival up to n — 1 is |(0|[/0„_i)p. To get the absolute 
probability for first detection in step n we have to multiply this with the survival 
probability With the normalization convention for we get the first 

arrival probability in the form |a„p, where a„ is the first arrival amplitude 

(10) a„ = (01 [/[/"- V), n>l. 

The total probability for events up to and including the n"^ step, i.e., detection at 
step fc < n or survival, thus adds up as 

n 

fe=l 

The return probability is therefore 

CO 

(11) i? = V|a„|2 = l- hm llC/'Vf. 

^ — ' n— f oo 

n=l 

Accordingly, we call the pair {U,4>) recurrent if i? = 1, and transient otherwise. 

This notion, which is the subject of our paper is not equivalent to SJK- recurrence, 
as the following example shows. 



RECURRENCE FOR DISCRETE TIME UNITARY EVOLUTIONS 



5 



Example 1. Let us look at an elementary example, which demonstrates the differ- 
ence between SJK-recurrence and the notion we introdnced. The Hilbert space will 
be 'H = C © ^^(Z), in which we wiU denote the basis vectors as |*) and \x), x gZ. 
The unitary is defined by U\*) = |*) and U\x) = \x + 1), and the initial state will 
be ^ = a\*) + /3\0) with |a|^ + |/3|^ = 1. Clearly, we have 

(12) (</.|C/"</)) = |a|2(l + M- 

To determine f/„0 we can solve the recursion for the expansion coefficients, which 
gives 

n 

(13) C/> = a|/3|2" I*) +I3\n)- I3\af J] IPf^"-'^ \n - k) 



Hence 

(14) s„ = ||C/>| 



This converges to s^o = 2|/3p/(l + |,5p), which is non-zero whenever ^ ^ 0. In 
other words, according to our definition, the only recurrent case in this example 
IS (j) = I*), which does not move at all. All other cases are transient. In contrast, 
p„ = |(</,|{7"</))|2 = |a|4 for all n > 1, so R^^^ = 1 unless a = 0. Hence SJK- 
recurrence holds with the only exception of the initial state which moves according 
the shift on ^^(Z). □ 

3. Spectral characterization of recurrence 

3.1. The mathematical setting. We have chosen our definition so that it only 

depends on the pair (U, 0) made up of a unitary operator on a Hilbert space "H and a 
unit vector (j) G H. It is clear from the definitions that only the subspace generated 
by the vectors C/"^ plays a role. We may thus ignore parts of the Hilbert space never 
explored by the evolution, and assume without loss that (j) is cyclic for U. Similarly, 
{U, (j)) are only needed up to unitary equivalence, i.e., the pair {VUV* , V(j)) with V 
any unitary operator yields the same return probabilities. One invariant of pairs 
{U, (f>) under unitary equivalence is the i/i-expectation of the spectral measure of U. 
This is the probability measure on the unit circle — {u G C : \u\ = 1}, defined 
by /u(rfu) = {(j)\E{du)(j)), where E is the spectral measure of U. Equivalently, we 
can characterize the system by the moments of n, i.e., by the Fourier coefiicients 

(15) Hn = jlJL{du)vJ' = {(I)\U''(I)), neZ. 

The measure /x is, in fact, a complete invariant: We may choose as a canonical form 
of the pair {U, (f)) the Hilbert space 'Hfj, = £'^{S^,n), the unitary C/^ of multiplication 
by the argument u, and the vector ^^(u) = 1. These are equivalent to the original 
setting by virtue of the operator V : H ^ "H^ defined by {VU"'<p){u) = u", for 
n e Z, which is unitary because {U"^(})\U"'(j)) = = Hn-m- 

Various functions have been introduced to characterize probability measures on 
the unit circle, all of which are analytic for \z\ < 1, and whose boundary values 
for 1 2; I — ^ 1 typically contain relevant information about the measure. We use the 
moment generating or Stieltjes function 



(16) ^(.) = f;^„.»= / 

n=0 



li{du) 
uz' 
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the Caratheodory function 

(17) F{z)^ ft,{du)'^ = 2ri{z)-l, 



1 F{z) - 1 _ 1 fi{z) - 1 
zF{z) + l ^ z f,(z) 



and the Schur Junction 

(18) f{z) = 

Here we have used the convention that, for an analytic function g, the analytic 
function with the conjugated Taylor coefhcients is denoted by i.e., g{z) = g{z). 
For \z\ < 1 we have 

(19) 5fteF(z) = / ^{du) = > 0, 

J \u-zY \\-zf{z)Y 

and hence |/(2:)| < 1. The Schur function thus maps the open unit disk to itself. 
The proof that, conversely every analytic function on the disc with ^eF > and 
F(0) = 1 (resp. |/| < 1) arises in this way from a probability measure /i is the 
subject of classic papers by Caratheodory, Schur, Hcrgiotz and many others. For a 
lively discussion of the history we recommend Simon's book [3D1 Ch. 1.3, Notes]. 

Both / and F have radial limits g(e**) = Muir^i- g{re^*) for almost all t so 
that we can consider that they are extended a.e. in the unit circle. The abso- 
lutely continuous part of /i is supported on the points z £ such that the radial 
limit satisfies 3ffeF(2) > 0, i.e., 1/(^)1 < 1. In fact, the density of the absolutely 



continuous part is given by (19). Therefore, the singular measures are character- 
ized by the fact that the related Schur function is "inner", i.e. \f{z)\ = 1 a.e. in 
|30j . The singular part of /i is concentrated on the points z E such that 
limr_>i- ^eF{rz) = oo, that is, zf{z) = 1. Concerning the pure point part of ji, 
the mass of any point z e 5^ is given by 

1 — r l — rl + rzf(rz) 
Li{]z\) = hm — - — F(rz) ^ hm — —. — -. 

" r^i- 2 ^ ' r^i- 2 l-rzf{rz) 



3.2. G ener ating functions. Let us now come back to the return problem as posed 
That is, we set U — (I — {(p) {(j)\)U , and ask for the survival probabilities 



Sect. 



2.3 



Sn — \\U"<j>\\'^ in the limit n — > oo, resp. the sum X^J^i |onP over the arrival 
amplitudes a„ = {(l)\UU^(j)) . 

To this end we compute a generating function for the vectors [/"(/i, starting from 
that of the operators C/", 

oo 

(20) G(z) = = (I-z[7)-\ |z|<l. 

This is essentially the resolvent of J7, and can be related to the corresponding ex- 
pression G{z) = (I — zU)~^ for the unitary evolution by the perturbation resolvent 
formula, which is especially efficient in this case, because the perturbation U — U 
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is a rank 1 operator. For completeness we repeat the short derivation here. We get 

G{z)-G{z) ^ (I - zU)-^ (^{I - zU) ~ {1 - zU)^! ~ zU)~^ 

= zG{zM{c^\UG{z) 

G{z){l + z\cj,){^\UG{z)) 
Giz)q^{l + z{<j)\UG{z)<j))) 
G{z)cp 



(21) 
(22) 



G{z) 
G{z)(t> 

G{z)(t> 



1 



z{cj,\UG{z)<j>)- 

Here, at ( |21[ ) the previous equation was apphed to (j). The function in the denomi- 
nator can be simphfied to 



(23) 



1 + z{(l)\UG{z)(j}) = L{du)(l + 



We now take the scalar product of (22) with U^(() to obtain a generating function 



for the arrival amplitudes, and use (23) also to simplify the numerator of (22). This 
yields 



a{z) 



(24) 
(25) 



n— 1 n— 

m(z) - 1 



n+l 



z{c^\UG{z)cf>) 



That is, the Schur function is essentially the generating function for the arrival 
amplitudes. This leads to an interesting analogy between the classical case and the 
quantum case discussed here: 





classical 


quantum 


return in n steps 
generating function 


probability 

Pn 

p see (|4| 


amplitude 
^ see (16 


1 


first return after n steps 
generating function 


probability 
q see ^ 


amplitude 
a„ 

a — zf see ( 


251 


renewal equation 


g = l-4 (§ 
p '-' 


1 i 

a = 1 - - ( 


24) 



In the classical case one takes the return probabilities Pn and uses the renewal equa- 
tion to get the first return probabilities . In our approach the renewal equation is 
used at the level of amplitudes, to pass from /x„ to the arrival amplitudes a„, which 
are then squared to give the arrival probabilities. This is in keeping with Feyn- 
man's elementary exposition of the "first principles of quantum mechanics" |12j. 
losely speaking: "In quantum theory add amplitudes and square at the end to 
get probabilities" . The continuous time limit of this idea is the root of his path 
integrals, and indeed there is an approach to continuous time measurements in 
quantum mechanics [21] based on the selection of suitable subsets of paths in the 
path integral. On the other hand, Feynman also explains that when you monitor 
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the process (his example is a measurement near the double slit), you should first 
square and then add probabilities. This would identify our approach as a "non- 
monitoring" one, in an apparent contradiction to our exposition in Sect. [2j It is 
debatable whether the principle "Add amplitudes where the classical theory would 
have you add probabilities" , which works so beautifully in Feynman's Quantum 
Mechanics 101, is really a fundamental law. In the end, it is Linear Algebra which 
can give us a unified understanding of the two columns in our table: Both columns 
are based on the powers of a matrix in which the first row has been replaced by 
zeros: in the quantum case by multiplying with the projection I — |'/>)((/)| and in the 
classical case by disregarding all jumps to 0. The renewal equation depends only 
on this, regardless of further structures such as unitarity and clcmcntwise positivity 
of transition operators. 

We remark that the above table suggests a reading of the SJK approach as "take 
the Pn from quantum mechanics, and define qn by the renewal equation" . This does 
lead to the SJK-recurrence criterion. However, in the paper '33| the quantities (?„ 
are avoided, as described above after The reason for this is not stated in the 
paper, but there is a good reason, which we discuss in Sect. |4.2| The g„ can be 
negative. 



3.3. The recurrence criterion. We can now utilize the identification (25) of the 
generating function a to turn the condition for recurrence, i.e. ^„ |anP 1, into 
a condition for the Schur function /. To this end we fix r < 1, and consider the 
series a(re**) = J2n ^"'^n 6**" as a Fourier series and use the Plancherel theorem: 

(26) ^r2"|a„|2^I^ C dt \f{re-'^')\\ 

Now recall that the Schur function is bounded by 1, so for the limit of this expression 
as r — > 1 to be equal to the maximum 1 we must have that |/(re**)| — !■ 1 for almost 
all t. In other words, / has to be an inner function. As we pointed out, this has a 
simple characterization directly in terms of the measure. 

Theorem 1. Consider a unitary operator U and an initial vector <j). Then the 
following are equivalent: 

(1) {U,4>) is recurrent. 

(2) The Schur function f is inner. 

(3) /I has no absolutely continuous part. 

Remark 1. Theorem [l] supplements known results about the interplay between the 
propagation properties of a state cf) and its spectral measure. For example a corollary 
of the Riemann-Lebesgue-lemma states that a state with absolutely continuous 
spectral measure will eventually leave any finite region, which corresponds to the 
statement of Theorem [T] that a state can only be recurrent and thereby returning 
to its initial starting point if its spectral measure contains no absolutely continuous 
part. In comparison to the famous RAGE theorem [22l[T5l|26] which characterizes 
the continuous and the pure point part of the spectral measure of a state by its 
dynamical properties our result distinguishes between its absolutely continuous and 
its singular part. 
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3.4. A spectral view of classical and SJK-recurrence. 



3.4.1. Classical. Here is brief account of the use of the "spectral method" to analyze 
classical random walks. Not all Markov chains with state space X can be analyzed 
in this fashion since the matrix P^^y can be very far from being symmetric. 

It turns out that the class of Markov chains that can be put in sclf-adjoint form 
has a nice description, namely they are the so-called reversible chains defined as 
those for which there is a vector with positive components tt^ such that 

Assuming the existence of this vector one can introduce an inner product between 
functions / and g defined on X by the rule 

xex 

One now checks easily that {Pf,g)-K = {f,P9)-K which says that the matrix P 
is now symmetric with respect to such inner product and the spectral theory for 
self-adjoint (bounded) operators can be used to great advantage. 

From the spectral measure E of our self-adjoint operator one gets a spectral 
measure m{d\) = {tl>,E{dX}i'j)T,. supported in [—1, 1], for any probability distribu- 
tion ipx- The probability of going back to the same probability distribution tjjx in 
n steps is given by the moments of the spectral measure, namely, 

m„ = jm{d\)X^ = {il^,P'^^)^, n = 0,1,2,... 

In particular, the probability p„ of returning to a given initial state after n 
steps coincides with the n-th moment of the spectral measure m for the distribution 
probability 5x.o concentrated at a; = 0. This identifies the generating function of 
the return probabilities j3„ as the Stieltjes function of the measure m, i.e., 



n=0 "^-1 



1 m{d\) 



Xz 

We finally have that the state is recurrent exactly when 

m(dX) 

Besides, using the renewal equation, we find that, in the recurrent case, the expected 
time for the first return to this state is 



^ d 



= lim — — 'ii^ = lim 



z=l 



>i 1-z z^i{l-z)p{z) to({1})' 



Therefore, the process returns to a state in a finite expected time exactly when its 
spectral measure has a mass at I, and the expected return time is the inverse of 
the corresponding mass. 

All this shows that, classically, the recurrence properties of a state only depend 
on the behavior of its spectral measure around the point I. 
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3.4.2. SJK. For the sake of comparison, we also give here a simple spectral char- 
acterization of SJK-recurrence, which is missing in |33| and was formulated in a 
different way in [5]. Since 



are just the Fourier coefficients of the measure, and /^_„ = /i„, the condition 



means, by the Plancherel Theorem, that fiidu) = p(u) du for a function p € 
{S^ , v) , where v denotes the uniform measure on the unit circle . Note that 
since is compact, this is more restrictive than just absolute continuity, i.e., the 
existence of an integrable density p. So we can summarize: A pair (U, (j)) Jails to 
he SJK-recurrent if and only if p is absolutely continuous with a square integrable 
probability density. 

3.5. Dependence on the starting point. One feature, which makes Polya's 
theory of recurrence so useful is that, if the walk is irreducible, recurrence is inde- 
pendent of the starting point, hence a property of the walk itself. Here "irreducible" 
means [18] that from any point every other can be eventually reached with non-zero 
probability, a property which is usually checked quite easily. 

In the quantum case the natural analogue of "the set of points reachable from 
starting point x" is "the cyclic subspace generated by the initial vector (/)", i.e., the 
span of all C/"0, which we will briefly denote by {U^cf))-^-^. Now it is clear that 
in the quantum case no non-trivial unitary operator can be "irreducible" , because 
the spectral decomposition will generate many subspaces which cannot be reached 
from each other. Nevertheless, some degree of independence can be stated, and 
follows immediately from Theorem [T] 

Corollary 1. Let 4>' G and {U,(j)) recurrent. Then {U,(f>') is recurrent. 

Note that this is not true for SJK-recurrence. As a counterexample we can take 
Example [ij or any other where p has both pure point and absolutely continuous 
components. This will be SJK-recurrent, but a suitable initial vector (f)' from the ab- 
solutely continuous subspace will be transient. A very similar argument shows that 
the above corollary is not true for transience. It is also not true for SJK-transience: 
this fails because even if (j) is absolutely continuous with square integrable density, 
{U^(p)-^-^ will also contain vectors whose density is not square integrable. 

We could also introduce a a notion of "hereditary transience" of {U, </>) meaning 
that all (j)' G (t/^ (/))-'-'- are transient. Clearly, this is equivalent to p being absolutely 
continuous. For hereditary transience the analogue of of the above corollary holds 
by definition. However, for the SJK version this is not an option (hereditary SJK- 
transience is vacuous). To complete the picture we might add that hereditary 
SJK-recurrence means that p has no absolutely continuous part, since that part 
would contain also vectors with square integrable density. Hence hereditary SJK- 
recurrence is exactly equivalent to recurrence in our sense. 

3.6. An inverse problem. A famous result in the classical case (37] states that any 
sequence g„ with g„ > 0, go = and qn = 1 can be the sequence of first return 
probabilities of a Markov chain. Here we discuss the analogous questions in the 



(27) 




(28) 




n— — oo 
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quantum case. It is clear from Bochner's theorem that the sequences /i„ — {(j)\U"' 
are exactly the positive definite ones with /ig = f ; the kernel kn^m = P^n~m has to 
be positive definite, where we set /x_„ = JI^ for n > 0. Equivalently, the sequence of 
Toeplitz determinants [30] is positive. Clearly, by expressing the /i„ in terms of the 
a„ this implies a hierarchy of conditions on the a„ . The first few look encouragingly 
simple: 

(29) l"|aip > 

(30) (l-laiH^-la^p > 

(31) (l-|ain3-(l-|ain(2|a2p + |a3n + |a2p > 23?e(a^aia3) 

These already suffice to answer some basic questions. First of all, we see that 
there is no quantum analog of the Yosida-Kakutani Theorem cited above: there are 



non-trival constraints on the sequence of first return probabilities. From (30), in 
particular, we see that the set of admissible sequences |a„p is not convex. Finally, 
the Toeplitz determinant constraints cannot be expressed in terms of the absolute 
values of the a„ alone. So even from the whole hierarchy of such inequalities it is 
not trivial to extract the precise constraints for the |a„|. 

4. Expected return time in the recurrent case 

4.1. It is an integer. Consider a recurrent pair ([/, 0). Its first arrival probabilities 
|a„p are then a normalized probability distribution on N. What are its features? 
Perhaps the first quantity to look at is the expected return time, i.e., 

oo 

(32) r = ^|a„|2n. 

It turns out that r can be computed from the measure very easily and turns out to 
be always an integer or infinite. This is surprising, because the a„ clearly depend 
continuously on the measure /i, so a small change in /i could be expected to also 
change r only a little. Our aim is therefore not just to prove this statement but to 
make it intuitive by giving r a topological interpretation. 

Theorem 2. Let {U, <j}) be a recurrent pair with spectral measure fx, Schur function 
f and expected return time t G M+ U {oo}. Then the following are equivalent 

(1) T < oo. 

(2) f is a rational function. 

(3) /I is equal to the sum of n < oo distinct point measures with non-zero 
weights. 

Moreover, t = n, and the polynomial degree of numerator and denominator of the 
Schur function is n — 1. 

We will prove the implications (3)^(2)^(1) emphasizing the narrative of ideas 
leading to the identification t = n. The slightly more technical implication (1)=>(3) 
is given in a separate, slightly more formal style. 



It is clear from (16) that Jl = P/Q is the quotient of two polynomials P,Q of 
degrees (n — 1), n, respectively. Then we can write / as R/P with R= [P — Q)/z. 
i? is a polynomial, because P(0)/Q(0) = 1, and hence P — Q has no constant term. 
Clearly, its degree is also n — 1. R and P have no common zeros, since these would 
be common zeros of P and Q, contradicting the observation that /i has n distinct 
poles on the unit circle. This proves (2). 
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Now observe that / is a rational function which, as an inner function, must have 



modulus one on the unit circle, and hence the generating function ( 25 ) of the arrival 
amplitudes 

oo 

(33) g(t) = e'7(e**) = ^a„e^"* 

n=l 

has modulus one for all real t. So g{t) winds around the origin an integer number 
w{g) of times as t goes from to 2tt. Locally we can write g{t) — exp(i7(t)), so the 
angular velocity is 

(34) da{t)^^=g(r)-dtg{t). 

ig(t) I 

Integrating this over one period t G [0, 27r], we get 2'Kw{g), so 



(35) w{g) 



— dt g{t) - dtg{t) = (na„) = r, 

•'^ n=0 



Here we used the Plancherel identity and the observation from ( 33 1 that the Fourier 
coefficients of {dtg)/i are na„. Hence we have a topological interpretation of the 
expected return time r, which already makes clear that it must be an integer. We 
remark that this connection was found in a totally different context by H. Brczis 
following a hint by I. M. Gelfand [7 . 

To determine this integer we use the Blaschke factorization, which decomposes 
a rational inner function into a product of the form 

n-l 

Oik - z 



(36) /(z) = /3[]^ 



k — l 

with \ak\ < 1 and |/3| = 1. The ak are exactly the zeros of / inside the unit circle, 
which by the reflection principle must be matched by poles at l/ok outside the 
circle. The upper limit of the product expresses our previous analysis of the degrees 
of the rational function /. Now the winding number of a product of several functions 
is clearly the sum of the winding numbers. Each of the Blaschke factors contributes 
1. Indeed, the winding number of a single factor clearly depends continuously on 
a, so must be constant, and we may set a = (i.e., f{z) — z) to evaluate it. 



Observing that the generating function ( 25 ) contains an additional factor z we get 
T = w{g) — n, proving (1). 

Before proving the remaining implication, let us see how a small additional point 
mass changes the Schur function and the arrival time distribution. An example is 
shown in Fig. [l] It is based on the sum of two point masses at m = 1 and u = i. We 
add a small point mass at u = — 1 by convex combination with a small weight e. 
This leads to an additional zero of / close to z = —1, and a full turn of the phase 
near that point. The interval in which this turn happens is of order e. Thus for 
fixed low n the contribution to the Fourier integral for a„ is negligible. For higher 
n it becomes visible, and the unit contribution to the expected return time comes 
from the large n Fourier coefficients. We will discuss this in more detail in Sect. [5] 



Proof of Theorem^ part (1)^(3): 

We prove this by contradiction, i.e., we show that for non-rational inner functions 
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Figure 1. Schur function for a sum of two and three point mea- 
sures. Bottom: unit circle with positions of the masses, represented 
according to weight. Middle: absolute value of /. Top: the phase 
of / on the unit circle. The ribbon is traced out by a vector indi- 
cating the argument of /. It has fixed length, is based on a unit 
circle (thick black) and orthogonal to it. Left panel: two point 
masses at 1, i with equal weight. Right: same mixed with a point 
mass near —1 with weight 0.01. The narrow spike near the light 
point mass does touch the zero plane: it is a zero of / near the 
circle with large derivative. 



/ we have r — oo. We combine the expression (26) with the observation that 
(37) T^-E^" 

fc=0 

increases monotonely to n as r — 1~. Hence, 

2„ ^ 1-^^^ i_^2|j(^g»t)|2 



(38) r(/)= hm$:^KP = hm - y 



1-r 



2 



(39) = l+lim^r.tl-|^(-^*)l^ 



r^l- 27r 7-7r 1 — 

We now consider a product of inner functions /i, /2, which is, of course, again an 
inner function. Then I/1/2I < |/i| and, consequently 

(40) t(/i/2) > r(/i). 

Suppose now that / can be represented as a convergent Blaschke product 



(41) /(^)=/5n 



,afc| 1 - Qffc z 

where the phase factors in the product have been introduced to make each factor 
positive at z = thus avoiding erratic behavior of conventional phases. Clearly, we 
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can take out the first iV — 1 factors as /i, for wliicli we know t(/i) — N. Hence 
t(/) > N for all N, so r(/) = oo. 

Not every inner function is a Blaschke product, but we can invoke a result of 
Frostman [T5| (see also [T?| Thm 6.4.]). It states that, given an inner function f{z), 
its Mobius transform 



(42) 



faiz) 



1 - afiz) 



is a Blaschke product for some a in the unit disk. Indeed, Frostman proved that 
this statement holds for almost all a (outside a set of null capacity) , but one a will 
be enough for us. Since f{z) ~ [a ~ fa{z))/{l — afa{z)), the two functions satisfy 
the estimate 



(43) 1-1/(^)1^ 



1 



|1 - afa{z)\'- 



a-\uz)n> 



1 



1 



(i-i/a(^)n. 



Hence from ( 39 ) 

(44) 



t(/)>1 + 



1 



1 



(t(/o)-1). 



Now fa is a Blaschke product. If it is finite, then / is rational and hence itself 
a finite Blaschke product, a case already covered. When it is an infinite product, 
then t(/q) = oo and hence t(/) = oo. I 



4.2. Expected arrival in the SJK scheme. Although absent in SJK papers, a 
definition of expected return time can be given for the SJK approach too. Let us 
begin with a seemingly natural approach which, however, turns out not to work. 
Recall Polya's Theory from Sect.[2j it allowed the determination of the probabilities 
Qn for first return in n steps from the probabilities p„ to be back at time n. Since 
the S JK-approach is based on applying Polya's criterion to the p„ as computed from 
quantum mechanics, we already have a candidate for (?„ and hence the expected 
return time via the renewal equation We get (compare also Sect. 3.4.1 ) 



(45) 



Polya 



dq{z) 



dz 



lim 



1 " g(^) 
1 - Z 



lim (1 — z)p{z) 

2— >-i 



If the measure consists of finite point masses the /i„ are given by exp(i0fcn) 
and we can further evaluate ( 45 1 to 



(46) 



Polya 



1 - 



Z 




We hasten too add that in the quantum case t^"^^^ has no operational meaning at 
all. What is worse, the "first return probabilities" q„ on which it is based need not 
be probabilities at all (see Example [2]). 

Therefore, in the spirit of , we build a notion of expected return time on the 
process described after ([s]): It is a Markov process with state space 0, 1, . . ., such 
that from a state n > one goes to with probability p„, and to n + 1 otherwise. 
The probability for first return to in n steps and the corresponding expected 
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return time are thus 



(47) 



,SJK 



„S,IK 



fe=l 



?SJK 



with (7g = and = Pi- As discussed above, this adds up to = R 

For a given pair (J7, (f) we thus get three sets of "first return probabilities" : the 
|a„p with the first return amphtudes a„ from (10), and two sets based on the 
probabihties p„ = namely q^^^ after (47) and after Sect. 



2.1 



The 

following example shows that all three notions can actually lead to different values 
for the same dynamical system. 

Example 2. In the Hilbert space consider as the unitary operator U the real 
rotation matrix by some angle 9 and choose the vector (j) — (1,0) as the initial state. 
The corresponding spectral measure consists of two mass points of equal weights 
1/2 at exp(±j0), which corresponds to the Schur function 

cos(^) — z 



(48) 



1 - cos(e')z 



Given the Schur function, we can calculate the three series of first return probabil- 



ities Or, 



SJK 



and qn- In Fig. [2] these three are shown together for 6 = 7r/4. In 
that particular case = (0, 1), so — {<P\U'^<i^) = 0, and p„ = 1, 1/2,0, 1/2, . . . 
repeated ad infinitum. Now this cannot be the return probability sequence of any 
Markov chain: If there is a non-zero probability to stay at site for one step, this 
probability must be non-zero for all n. So p2 ~ Q is impossible classically. In the 
jargon of quantum information theory the unitary we have chosen is called a square 
root of the NOT gate [6 , since swaps the two bit values, a feat which cannot 
be realized by two equal classical probabilistic steps. Indeed, q2 — —1/4, and the 
sum of all the negative qn is around —0.77. It is also clear from Fig. [2]that there is 




10'- 
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20 
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Figure 2. The figure shows the first return probabilities (left 
panel) and the accumulated first return probabilities (right panel) 
for 9 = 7r/4. The three curves correspond to the first return prob- 
abilities given by |a„p (thick), g^''^ (long dashed) and qn (short 
dashed). 



no general inequality connecting the first returns in the SJK approach and in our 
approach: which of the two gives higher probability for the event "arrival before f 
may depend on t. 
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R 

2- — 



Figure 3. Dependence of the expected return times r = r^"'^^ 
(thick) and t^-^^ (dashed, first 400 summands) on the rotation 
angle 9 (for 6 e [— tt, 0] the graph shows the same behavior). 



To see that this is also true for the expected return times we evaluate the expres- 
sions for T, r^°'^'^ and r^''^ for this example. Since the overlap oi (f> = (1,0) with 
the two eigenvectors of U is independent of the rotation angle 6, Theorem [2] tells us 
that r equals the number of distinct eigenvalues of U. This amounts to r = 2 for 
all 9 except the singular cases and ±n, where U has two equal eigenvalues and 
therefore r is 1. Since both mass points have equal weights 1/2 equation (|46]) for 
-,-Poiya gives the same result. 

which, since p„ = cos^(n6') holds, 
amounts to 

n 

n>l k=l 

Fig.jsjshows the two curves of the three expected return times t, t^"'^'* and t^-^^ 
depending on the rotation angle. By varying 9, we find that r < r^"'^ as well as 
^SJK ^ ^ realized. In addition by choosing an initial state that does not 

have equal overlap with both eigenvectors also t^"'^'' < t^-^^ can be realized. □ 

We have seen that a single point measure contribution in /i makes a process 
SJK-recurrent. The following proposition makes an even stronger statement. 

Proposition 1. Let ^ be a measure on the unit circle with at least one atom. 
Then r^-'^ < oo. 



To determine r we have to evaluate (47 



Proof: Let s„ = Ofe=i(l ~ Pk) be the survival probabilities in the SJK-approach. 
Then qji -f Sn — Sn—1 and, by partial summation, 

m rn — 1 m m— 1 m— 1 

^nqn= + l)sn -^nsn= ^ s„ - msm < ^ s„- 

n— 1 n— n— 1 n— n— 

Using the inequality > 1 + x, a: G M, we find 

n 

^SJK < 1 + ^ [](1 - < 1 + ^ e-5:^^i \^^-\\ 

n>lA;=l n>l 



RECURRENCE FOR DISCRETE TIME UNITARY EVOLUTIONS 



17 



Wiener's theorem on the unit circle states that [3T] 

hm = hm — — = Hpp{S ), 

n— i-oo n n— >-oo Zn + i 

where fipp is the pure point part of the spectral measure. Hence, the presence of 
mass points implies that 

^^-1 > M > 0, n > 1, 
n 

for some constant M. As a consequence, 



n>l 



5. The finite case 



In this section we consider a finite dimensional system, so the measure is a 
sum of point measures at positions Ui with weight rui, i = 1, . . . ,n. Thus 



(49) ni^) 



^ 1 



and the Schur function / is a finite Blaschke product ( 36 ) . Our plan in this section 
is to give a more detailed account of the changes in the return probability when n 
changes and hence the expected return time t — n jumps. In particular, we will 
show that at these points the variance of the first return diverges. We begin with 
a formula for the variance in terms of the zeros Zj of /. 

Proposition 2. Let f{z) = X^fc'^-^*^^^ rational Schur function with zeros 

zi, . . . , Zn^i. Then 

(50) v. = Y. - (E - E 

k k j,l ^ 



Proof : We took the connection between / and the a„ from our generating function 
relation (25 1. For the sake of this proof, however, it is easier to omit the complex 
conjugate, which is anyhow irrelevant in (50 1, and to shift the sequence a„ so that 
the o„ are really the Taylor coefficients of /. Note that J^k '^kfcP = n — 1 is then 
the number of zeros. 

From the Blaschke product (|36| we get 

/' ^ x^/ -1 -zl 

^\z,-z 1 



(51) 



/ 





z7z 



Since |/| = 1 on the unit circle, we can express the second moment as 



fe 



2tt 



dfje^') 

de 



^""dO 
2ti 



f 
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where the circle indicates complex integration around the unit circle. Writing the 
absolute value as XX we combine the factor X with the factor 1/z. Then, using 
jzp = 1, all denominators are linear in z. This gives 

dz ( 1 Z£ \ / 1 zl 



(53) 5:fc2Kp = 5: 



^ 27ri \1 — zg z z — zg / \ z — Zj 1 — z,- _ 

k j,i 

The integral can be solved by residue inspection. The first term in the first paren- 
thesis and the second term in the second parenthesis are analytic in the disc. The 
product of the other two terms gives a zero integral when Zj — zg because the pole 
is then of second order, and also gives zero otherwise because it is proportional to 
l/(z — Zj) — l/{z — ze). Combining the remaining terms gives 

(54) Y.'^^iakf - E(r^ + - E 

k U "^""^ l-Z^ZtJ ^^\-ZtZj 

Now in order to get the variance, we have to subtract (n— 1)^ from this expression, 
or 1 from each term. This gives the expression stated in the proposition. ■ 



The expression (501 is positive, because we can write it as 

2 

11 — J. 

(55) K = 2'^ 



k=l 



This can vanish only if J^k '^k ~ ^ ^'^^ ^- since by Newton's relations 

[25] the coefficients of the polynomial Yiki^ ~ ^fc)' elementary symmetric 

functions of (zi, . . . , z„„i), can be expressed in terms of the sums of powers, the 
equation for the Zk can only be = 0. In other words, Vr = implies that 

/(z) = /3z"~^, = 1, corresponding to n equal point masses at the n-th roots of 
(3 e S^. This should not surprise us: It corresponds precisely to a "clock" unitary, 
which cyclically rotates n basis states, of which one is the initial state. Clearly, this 
returns to the origin after exactly n steps. For fixed n, small variance still implies 
that all Zk are small. 

In the opposite direction, suppose that one of the zeros, say zi, approaches the 



unit circle, while the others are fixed. Then the term with ^ = j = 1 in (50) 
diverges, while all other terms stay finite. Hence Vr oo. Characteristically this 
happens when the number n of point measures, and hence the number of zeros 
changes its value. Indeed this could happen either by fusion, when two points 
converge to each other (e.g., u„ — ^ Wn-i), or by a shrinking, when the weight of 
one point goes to zero (e.g., to„ — 0). In either case we get a weakly convergent 
sequence of measures for which n changes in the limit. But by weak convergence the 
sequence of Schur functions converges uniformly on every disk of radius (l — e). This 
implies that the number of zeros of the Schur function in any such disc converges. 
In other words, the change in the number of zeros can only happen by one zero 
converging to the unit circle, which entails — oo, as we have seen. When the 
change in the measure is localized near one point, as for fusion and shrinking, then 
the uniform convergence argument also holds for a neighborhood of other parts of 
the unit circle, and we can even say that the vanishing zero has to converge to the 
point with shrinking weight (as shown in Fig. [ij or to the fused point. 

One source of intuition for such connections between the point measures and the 
Schur zeros is the following electrostatic interpretation. Although it is stated for 
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finite sums of mass points, which is the case under study in this section, it can be 
easily generalized to infinitely supported measures too. 

Lemma 1. Consider a measure fi on the unit circle, which is a finite sum of point 

measures, so 'il{z) — X)r=i™i/(1 ~ '^i^) with rui > 0, — 1; '^'^'^ I""*! = 1- 

Then the zeros Zi, . . . , 2;„„i of the Schur function are precisely the stationary points 
of the two-dimensional electrostatic potential 

n 

(56) y(z) = ^TOaog 

i=l 

In particular, all Zk lie in the convex hull of the Ui . 

The mi are the point charges of this electrostatic system located at the Ui. So 
each of the point masses generates a repulsive force field. If one of the charges 
is small, there will be a stationary point close to this charge at Ui, as seen in Fig. [T] 
The convex hull of the Ui is indicated in red in Fig. [l] 

Proof : Clearly, by ([Ts]) the zeros of / are those of 

-U{z)-l \ = -> m, ^-1 = > T — = > . 

Z\ / Z ^ — ' V 1 — ULiZ I — ' 1 — UiZ — ' Ui — z 

It is easily verified that this is equal to 

dV{x + iy) . dV{x + iy) 
dx dy 

so at the zeros of / the gradient of V vanishes. The statement about the convex 
hull may seem obvious from this. More formally, suppose to the contrary that the 
above expression vanishes for some z outside the convex hull of the Ui . Then there 
must be a linear functional z M- 5Re(Az) and a constant c such that ?Sie{Xui) < c for 
all i, but die{Xz) > c. But then 

which is a contradiction. ■ 

We note that the arguments which we gave for the connection between low 
variance of first return times and concentration of Schur zeros near the origin is 
valid only for fixed n. i.e., explicit bounds, which could be derived from these 
arguments will depend on n. This is shown by the following example, in which all 
zeros (an increasing number) go the unit circle, and the variance goes to zero. We 
choose the zeros to be Zj = Xe^'^'^^^ ("-i) , for fc = 1, . . . , n — 1. Then in the sum ( 55 ) 



the sum in the absolute value is non-zero only when s is a multiple of {n — 1), in 
which case it is equal to (n — 1)A''. Summing over these multiples and noting that 
s = is excluded also, we find 

(57) V. ^ ______ . 

Hence if we choose logjAp = — 3(n — l)^^log(n — 1) we have |A| — > 1 and x 
(n- l)-i ^ 0. 
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6. Examples with absolutely continuous measures 



In this section we deal with transient pairs (U,4>), i.e., those whose spectral 
measure fi has an absolutely continuous part. We focus on ways to get the first 
return probabilities |a„p or at least their sum, the total return probability 



A typical class of such models are so-called quantum walks, i.e., systems with a 
discrete spatial degree of freedom, which in each step can change only by a finite 
distance. There is a new interest in quantum walks coming from recent experimental 
realizations [TU]. These systems were also the motivation for [33", as well as for us. 
This section is therefore split between two approaches to quantum walks, which 
have their characteristic strengths, but also a large overlap. 

The first method, introduced by two of us (F.A.G. and L.V.) with coauthors 
[51 [3] builds on a method turning any pair {U,(f>) into a quantum walk. This walk 
most naturally lives on Z+ ~ {0, 1,2,...}, although it can be extended to deal with 
Z. It almost always has position dependent steps. When the dynamics is periodic 
from some site onwards, one can determine the Schur function by a fixed point 
equation, which typically gives absolutely continuous spectrum plus some discrete 
eigenvalues related to the initial segment. 

The second method is designed for strictly translation invariant systems in any 
lattice dimension. It has been introduced in [TB] and was used by two of us (A.H.W. 
and R.F.W.) with coauthors [1] for getting the asymptotic position distribution 
even in the presence of time dependent noise. The core of the method is the 
joint diagonalization of U and the translations, which amounts to diagonalizing a 
momentum dependent matrix. Finite range spatial variations in the dynamics can 
be dealt with via perturbation theory , which again produces additional discrete 
eigenvalues analogous to the introduction of impurities into a periodic potential. 

Both methods apply to so-called "coined walks" in one dimension, which live on 
the Hilbert space ^^(Z) (g) C^, whose basis vectors we write as \x,a), x G Z, and 
a E {t, i}- The time step unitary is factorized into "shift" and "coin" U = SC with 
S\x, t) — |a;-|- 1, t) and S\x,l.) = \x— 1,1), and a coin operation, which acts at each 
site X separately as a 2 x 2 unitary matrix Cx on |a;,^). The translation invariant 
case corresponds to a site independent coin Cx = Cq, and will be treated in both 
approaches. In contrast to the classical result of Polya, who found recurrence for 
one dimensional unbiased nearest neighbor random walks the coined quantum walks 
are always transient. Intuitively, this may be related to the much slower spreading 
(~ rather than ^ t) of the classical walks: they just spend more time close to 



6.1. Return probability from orthogonal polynomials. How can we turn any 
pair {U, 4>) into a quantum walk? The basic idea is to build a basis of the Hilbert 
space from the Gram-Schmidt orthogonalization of the sequence 4>, U~^(j), U(j), 
U~'^4', . . ., which is equivalent to building the //-orthogonal Laurent polynomi- 
als on the unit circle. This leads to a canonical five-diagonal matrix representation 




oo 



home. 
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of [/, a so-called CMV matrix [TOl |30l |32l |35] 



C = 



1 % 


Po 














Po7i 


-7o7i 


pi 72 










PQPI 


-7oPi 


-7i72 


-71 P2 














P273 


-7273 


P374 


P3P4 








P2P3 


-72 P3 


-7374 


-73P4 














P475 


-7475 














PApT, 


-74P5 















l7fe| < 1, 



Pk 



lk\ 



Then the pair (C,eo), with eg = (1,0,0,...)* the first canonical basis vector is 
unitarily equivalent to (U,4>). 

The coefficients 7^ are called the Schur parameters of / or the Verblunsky coef- 
ficients of the measure ^ [Ml [3D] . They yield a continued fraction expansion of the 
Schur function / through the so called Schur algorithm [25] 



(58) 



h{z)^f; fk+iiz) ^ 



1 fkjz) - 7fc 

2 1 - Ikfkiz) ' 



lk = .fk{0). 



Each fk is a Schur function in its own right, which has Schur parameters 7^, 7^+1, 
7fc+2, • ■ • ■ According to a remarkable result due to S. Khrushchev [2l] it is closely 
related to the pair {C,ek) starting at the basis vector e^, which has the Schur 
function fkBk, where is a finite Blaschke product. Hence is equal to the 

total return probability starting from Cfc. 

Coined walks on Z-|_ are a special case of this general construction, namely the 
case that the odd Schur parameters 72^+1 vanish. Thus we get the CMV matrix 

\ 



(59) 



/ 

Indeed, if we relabel our basis as e2x = \x,t) and 62^+1 = we get precisely 

the coined quantum walk on the half line with reflecting boundary condition and 
the coins 
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(60) 



P2x 

l2x 



-72a 
P2x 



\l2x\ < 1, P2x = - l722;|" 



Conversely, given a sequence of coins, we can immediately build the CMV repre- 
sentation C of the corresponding walk. We will assume from now on that |72£c| < 1, 
which makes the walk "irreducible", i.e., there is no finite interval in x which is 
decoupled from its complement. 

Because the odd Schur parameter vanish the iteration (581 gives /2£!;(z) = z^^ f2x-i 
We conclude that the states |a; — 1, |) and |x, t) of an irreducible two-state quantum 
walk on Z_|_ have the same return probability, and this return probability does not 



depend on the coins Cj,, y < x. 
in a separate publication. 



This is a remarkable feature, which we will explore 
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6.1.1. ID Quantum walk with a constant coin on Z+. As an example we will con- 
sider a quantum walk with a constant coin of the form (60) with 722: = 7 constant. 
The Schur parameters 7,0,7,0,7,0,... and hence the sequence of Schur iterates 
fk is hence periodic with period L = 2 |5J I3 ■ Therefore it satisfies a fixed point 
equation for a composition of L steps of Schur iteration. In our case / = /2 is just 
a quadratic equation, which has the solution 



(61) /(^) = 



27^2 

The branch of the square root is that one converging to 1 when z — )■ 0. 

The measure related to /, i.e. the spectral measure of the cycHc vector |0,t), 
has a weight and a mass point which is absent only when the coin is symmetric 
[5]. Hence, all the states are transient in case of a symmetric coin, otherwise a 
one-dimensional transient subspace appears, namely, the eigenspace of the simple 
eigenvalue given by the mass point. 

The return probability for the state |0,t) is 

-it 

./n 27r 7 



sign(cosi) \ \^\'^ — s,\v? t if |sini| < I7I, 



—sign (sini) iJ siv? t — |7|2 if | sini| > I7 



where 

h{t) = 
The final result is 

(62) = -1^ {pI^I + (1 _ 2p2)^} , p= 77 = arcsin|7|, 

which actually holds for any basis state due to Khrushchev's result and the fact 
that f[z) — f2x{z) = f2x-i(,z). The return probability does not depend on the 
phase of 7, and it is an increasing function of |7|. 

Using the generating function of the Legendre polynomials P„, 

^ = y p„(x)z", 

VI - 2xz + z2 'X'. 

n>0 

we can also find the power expansion of /, obtaining for |0, t) the arrival amplitudes 
a„, where ai = 7 and 

n Pn-l\X) — xPn\x) -, m |2 \ i 

a2n = 0, a2„+i = — — — , a; = l-27 , n>l. 

2j[n + 1) 

Bearing in mind that P„(cos(0)) ~ ■^2/7rnsin(6') cos[(n + 1/2)9 — 7r/4] we find that 
the first return probabilities |a„|2 decay as n^'^, 

|a2„+iP ~ ;r^{l - sin[(4n -I- 2)r,]} 
27r|7| 

6.1.2. ID Quantum walk with a constant coin on Z. Following the approach in 
[Eli, we choose the basis = {|0,t), |-l,i), |-l,t), |0,;), |l,t), |-2,i), 

|— 2,t), |l,i), . . . } to mimic the half infinite situation. This leads to a CMV repre- 
sentation C of the unitary step operator similar to ( |59[ ) with vanishing odd Schur 
parameters, but matrix valued, i.e., 7^ and pk are now antidiagonal respectively di- 
agonal 2x2 matrices. The Caratheodory function and the Schur function likewise 
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become matrix valued. The final result is that, for any basis state, the first return 
probabilities |a„p decay as like in the previous case, 

2p 



a,in-\ = 0, |a2„| 



{1 + sin(4n7])}n ^, p= - Yl\^ , 



T] = arcsm |7| 



and the total return probability is given in terms of the Schur function / in (61), 

2 



(63) 



R^\\f\ 



2||2 



{(l + 2p2)p|^| + (l_ V)^}. 



We have a similar behavior as in the case of Z+ but, for the same coin, the return 
probability is slightly smaller for Z than for Z+, as Figure |4] shows. 




Figure 4. The return probabihty R for a basic state |A:,t), 
of a quantum walk with a constant coin C — {cij ) as a function of 
I7I = 1^21 1- The upper curve corresponds to a constant coin on the 
non-negative integers, the lower one on the integers. 



We can summarize stating that the recurrence properties for a constant coin 
on Z or Z_|_ are similar, but the recurrence is slightly more prominent when the 
translation invariance is broken by the boundary conditions in Z+. 

6.2. Return times by the Fourier method. This method was developed for 
translationally invariant quantum walks [161 0]. In this sense it is much more 
special than the CGMV method described above, but has the advantage that it 
works just as well in higher lattice dimension, and with an arbitrary initial state, 
without the need to recompute the Verblunsky coefficients. Moreover, the Fourier 
method has been extended to certain decoherent walks and interacting two-particle 
systems [llllll]. 

The Hilbert space is now ^^(Z'*) K., where s G N is the dimension of the under- 
lying spatial lattice, and JC is the space of internal degrees of freedom at each site 
X. We only assume that U is unitary on this Hilbert space, and commutes with the 
lattice translations. Often it is made part of the definition of a quantum walk that 
it has strictly finite propagation in each step, but we will not need this condition 
here. Therefore, we can jointly diagonalize U with the translations by a Fourier 
transform in the spatial variable, which means that the Hilbert space vectors are 
naturally considered as /C-valued functions of the Fourier variable "momentum" 
P = G [^^j""]^- In this representation U becomes multiplication by 
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a p-dependent unitary operator on /C, denoted by U{p). The assumption of fi- 
nite propagation speed would imply that the matrix elements of U{p) are Laurent 
polynomials in the variables exp(ipfc). 

Then for any initial vector (jj we get the Stieltjes function ([l6| 



(64) 



When the starting point is the origin, 0(p) = G /C is independent of p, and we 
can evaluate the integral as /I(z) = {(f>\M{z)\(f)) , with the Stieltjes operator 



m(z) = (01(1 - zt/)-V> = / dp (0(p)|(l - zU{p))-'\^ip)) 



(65) 



M{z) = dp {I- zU{p))- 



For a coined one-dimensional quantum walk with the same coin as in Sect. |6.1.1| 
the walk operator takes the form 



(66) 



U{p) = S{p) ■ C = 



1 P 



The integral ( 65 ) can be evaluated by the residue theorem and gives 



(67) 



M{z) 



1 / 1 - + g{z) 
2^) \ l{l + z^-g{z)) 



-lil + z^~g{z)y 



l-z' 



with g{x) = \/ — z^)^ + 4|7pz2. From this we can immediately write down the 
Schur function and compute the total return probability R as the 2-norm. In 
agreement with Sect. 6.1.2 the cases ^ = | f) and (fi — \ I) give the same result (631, 
because the diagonal elements of M{z) are equal. However, for starting/absorbing 
states (j) which are superpositions of these two we get other Schur functions, because 
M{z) is not a multiple of the identity. 



0.41 r 




Figure 5. ^-dependence of the return probability i? of a quantum 



walk with 7 — l/v2 (see (66 1) of the initial state (j) = (l:exp(z0)). 



Indeed, this can lead to different R. This can be seen from fig. [5] which shows 
the dependence of R on the relative phase of the initial state </> — (l,exp(i6')) for 
the quantum walk with 7 — l/\/2 according to (66). 
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6.2.1. Local perturbations. In the CGMV approach a local perturbation is easily 
implemented by letting the Schur coefficients be constant only from some index 
onwards. From there on the fk are determined by the fixed point equation, and 



/ = /o is obtained by inverse Schur iteration ( 58 ) . 

In the Fourier approach a local perturbation can be turned into a perturbation 
only at the origin, by redefining a block of cells as a single one. The resolvent 
of the perturbed unitary operator can be evaluated with the resolvent formula for 
perturbations. Since the "unperturbed" resolvent of the above operator U can 
be computed separately for each U (p) this leads to a finite dimensional, albeit p- 
dependent matrix computation. Of course, as in the CGMV case this may quickly 
lead to rather unwieldy expressions. However, at least for numerical work one gets 
a systematic method. 

7. Outlook 

7.1. Examples with singular measures. Obviously, there is a wide range of 
examples which one would like to treat under this heading. In the quantum walk 
field, disordered walks, i.e., operators with space dependent i.i.d. random coins come 
to mind. In this case one has Anderson localization [3|, i.e., dense point spectrum. 
In almost periodic examples, such as walks in an irrational magnetic field, one 
expects purely singular spectrum. For irrationals which are well approximated by 
rationals, e.g., Liouville numbers, there may be late returns grouped in a hierarchy 
of time scales. A further family of examples is given by measures with some self- 
similarity such as a famous example by Riesz [571133 HZ] -For the time being we have 
no general results about return times in these examples, but it seems to be a rich 
field for investigations linking spectral properties with the dynamics of propagation, 
which we intend to explore further. 

7.2. Generalizations. We have focused completely on the case where absorption 
happens only in one pure state, which is also the initial one. The main reason 
for this was to have a very direct link to spectral theory, and a dependence only 
on the pair (t/, </)) and nothing else. However, from the point of view of physical 
applications some generalizations are desirable, and can be pursued very much 
along the lines set out in this paper. Indeed there is no reason why in Sect. |2.3| 
the initial and the absorbing state should be the same, so one can immediately 
generalize to "first arrival" rather than "first return" . It is often natural to consider 
a more than one dimensional space to be absorbing, e.g., when we want to discuss 
the return of a coined walk to the origin, regardless of the internal state is. In 
that case one might also want to look at spin resolved arrival events and, quite 
generally, at an array of counters. This has been studied in the continuous time 
case [51 which is a further generalization of interest. When an appropriate 
weakening of the absorption process is introduced to avoid the Zeno effect [11] 
in the continuous time limit one finds a picture analogous to the discrete time 
case: The absorbing perturbation then modifies a one-parameter unitary group to 
a semigroup of contractions. The difference of their generators, sometimes called 
an optical potential, may be unbounded or even fail to be a closable operator, as in 
the case of an absorbing point "potential" —iXS{x) perturbing the Laplacian (free 
particle) . The interplay between spectrum of the Hamiltonian and arrival time in 
the continuous case contains a recently established uncertainty relation [20] . but 
much is left to explore. 
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